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Инструментальные средства оперативного анализа 
местности для построения систем 
принятия решений при чрезвычайных ситуациях 


В статье представлены разработанные инструментальные средства оперативного анализа местности на 
основе цифровых карт и снимков для решения задач по моделированию, анализу и принятию решений при 
чрезвычайных ситуациях. Инструментальные средства поддерживают решение таких задач как анализ 
свойств местности, расчет кратчайших маршрутов с учетом заданных ограничений, поиск ближайших 
объектов с заданными характеристиками, определение пространственных отношений объектов местности, 
расчет зон затопления, моделирование переноса радиоактивных частиц, моделирование распространения 
пожаров и ряд др. Приводятся примеры использования. 


Введение 


Для оперативного и качественного принятия решения и выполнения мероприятий 
по ликвидации последствий чрезвычайных ситуаций необходимо оперировать с инфор- 
мацией, имеющей следующие основные свойства: 

— обеспечение высокой точности локализации места ЧС; 

— обеспечение оперативного получения актуальной информации о местности, 
на которой произошло ЧС; 

— наличие достаточных характеристик свойств местности, обеспечивающих ка- 
чественный анализ возникающей ситуации. 

Такими свойствами в полной мере обладают цифровые карты местности (ЦКМ) 
и данные дистанционного зондирования Земли (ДЗЗ). 

В ОИПИ НАН Беларуси разработаны инструментальные средства оперативного 
анализа местности на основе цифровых карт и снимков для решения задач по моделиро- 
ванию, анализу и принятию решений при чрезвычайных ситуациях. В качестве базы для 
построения программного обеспечения положены модели и методы оперативного ана- 
лиза местности для поддержки принятия решений на основе цифровых карт и снимков, 
изложенных в [1]. Инструментальные средства включают решение таких задач, как 
анализ свойств местности, расчет кратчайших маршрутов с учетом заданных ограниче- 
ний, поиск ближайших объектов с заданными характеристиками, определение прост- 
ранственных отношений объектов местности, расчет зон затопления, моделирование 
переноса радиоактивных частиц, моделирование распространения пожаров и ряд других. 
В данной статье приведены методы решения задач нахождения оптимального маршрута 
движения транспортного средства с учетом ограничений на проходимость и присутствие 
зон поражения и моделирования переноса радиоактивных и загрязняющих частиц, ко- 
торые могут поступать в атмосферу в результате взрыва, пожара, выброса примесей при 
аварии на опасном объекте, в результате военных действий или атак террористов, и 
описываются программные средства, реализующие эти методы. 
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Поиск кратчайшего расстояния 


Программа осуществляет расчет наиболее выгодного маршрута перемещения транс- 
портного средства по наличной дорожной сети. Критерии выгоды могут быть по времени, по 
пройденному расстоянию, по условиям проезда. При прокладке маршрута должны учи- 
тываться проходимость отдельных участков дорожной сети, крутизна подъемов и спусков, 
наличие мостов, проездов и других ситуаций, ограничивающих проезд. При необходимости 
можно задавать режим объезда населенных пунктов и заданных зон на местности. 

Входные данные: 

а) карта местности с тематическими слоями: дорожная сеть и сооружения; рельеф 
местности; растительность; 

6) населенные пункты с характеристиками численности населения; 

в) полный граф дорожной сети для этой местности; 

г) прореженный (упрощенный) граф дорожной сети, полученный из полного 
графа, путем исключения из него неключевых узлов; 

д) характеристики транспортного средства, для которого прокладывается маршрут; 

е) характеристики дорожной сети; 

В общем виде решение сетевой задачи состоит из двух этапов: 

— построение графа сети; 

— нахождение оптимального маршрута на графе. 

Построение графа сети 

Для построения графа дорожной сети используется следующий алгоритм: 

1. Из карты выделяется слой дорог и формируется массив дорог. При этом 
можно заранее ограничить объем расчетов, исключив из массива. 

2. Создается пустой граф дорог. В качестве вершин графа выступают опорные 
точки дорог и перекрестки дорог. В качестве ребер графа выступают прямые отрезки 
дорог, соединяющие опорные точки. 

3. Формируется цикл по массиву дорог. 

4. Формируется вложенный цикл по опорным точкам дороги. 

5. Производится проверка на независимость новой вершины графа. При этом 
могут возникнуть следующие ситуации: 

1) новая вершина не имеет в пределах некоторой окрестности уже присоединен- 
ных вершин графа или ребер и новое ребро ни с чем не пересекается и не примыкает. 

В этом случае к графу присоединяется новая вершина, которая присоединяется 
к предыдущей вершине ребром. 

2) новая вершина находится в некоторой окрестности или совпадает с одним из 
узлов графа. 

В этом случае новая вершина в графе не создается, а новое ребро присоединяется 
к уже существующей вершине. 

3) в некоторой окрестности новой вершины имеется существующее ребро. 

В этом случае к графу добавляется новая вершина и новое ребро, а существующее 
ребро удаляется из графа, вместо него добавляется еще 2 новых ребра, которые соеди- 
няют новую вершину и вершины, которые соединялись существующим ребром; 

4) в некоторой окрестности нового ребра существует вершина графа. 

В этом случае к графу добавляется новая вершина и два ребра — от предыдущей 
точки к существующей вершике, и от существующей вершины - к новой. 

5) новое ребро и существующее пересекаются. 

В этом случае к графу добавляется новая вершина, в точке пересечения ребер со- 
здается еще одна вершина, которая добавляется к графу, существующее и новое ребро 
удаляются из графа, а вместо них создаются 4 новых ребра, которые соединяют 4 вер- 
шины со сформированной точкой пересечения. 
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После этого переходим к следующей точке дороги, или новой дороге. 

В случае если новая вершина является первой точкой в рассматриваемой дороге, 
то к ней применяется только ситуации 1,2 и 3. 

Оптимальная величина окрестности для проверки независимости вершин и ребер 
равна 5 — 10 дискретов оцифровки карты. 

После создания полного графа связности дорожной сети по этому графу фор- 
мируется, если необходимо, тематическая информация: наличие мостов, туннелей, над- 
и подземных переходов, пересекающих линий электропередач и путепроводов, вели- 
чина уклонов, тип дороги и прочее. 

После сбора всей необходимой информации создается упрощенный граф дорож- 
ной сети. Упрощение заключается в следующем: 

а) в упрощенный граф заносятся только вершины, которые имеют одно (концевые 
участки) или более двух (перекрестки) ребер; 

6) создаются ребра упрощенного графа, каждое из которых является суммой ребер 
соединяющих через проходные вершины полного графа; 

в) на ребра упрощенного графа переносится тематическая информация ребер 
полного графа. 

Трудоемкость данного алгоритма при создании графа дорог составляет О(п х 10°(п)), 
где п — количество точек дорог. Поскольку количество точек в реально существующей 
карте может составлять порядка 50 — 100 тысяч, то рекомендуется созданный граф со- 
хранять в виде отдельного картографического слоя. 

Нахождение оптимального маршрута в графе 

Следующим этапом после создания графа дорог является нахождение кратчай- 
шего пути между двумя вершинами в графе со взвешенными положительными ребрами. 

Решение задачи определения оптимального маршрута выполняется методами типа 
Дейкстры, Беллмана-Форда. 

Для построения графа дорожной сети необходимо задать следующие параметры: 

— величина погрешности — максимальное расстояние между отдельными объектами 
дорожной сети, при котором объекты склеиваются. Задается в дискретах оцифровки карты; 

— классификационные коды дорог, включаемые в состав графа, — список классифи- 
кационных кодов, которые задают типы дорог для включения в граф. 

Процесс создания графа в зависимости от размера карты и выбранного массива 
классификационных кодов может занимать от секунд до часов. Поэтому имеет смысл 
хранить в виде отдельного картографического слоя. 

Для поиска кратчайшего пути пользователю необходимо загрузить карту и граф 
дорожной сети, после чего задать начальную и конечную точку пути и вызвать функцию 
расчета (рис. 1). 


Линейное 


препятствие 


препятствие 


Рисунок 1 — Поиск оптимального маршрута по дорожной сети с оперативным 
указанием препятствий 
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Моделирование переноса радиоактивных 
и загрязняющих частиц 


Цель задачи — разработка функций для прогноза ветрового переноса сильно дейст- 
вующих ядовитых веществ, аварийно химически опасных веществ, радиоактивных и пас- 
сивных загрязнений при чрезвычайных ситуациях. Загрязнения могут поступать в атмо- 
сферу в результате взрыва, пожара, выброса примесей при аварии на опасном объекте, в 
результате военных действий или атак террористов. Выбрасываться могут радиоактивные, 
химически опасные, твердые и жидкие вещества в виде дыма и пыли. 

Источники выбросов могут быть стационарными и подвижными, точечными, ли- 
нейными и площадными. Линейные и площадные источники трактуются как совокуп- 
ности точечных источников. 

Расчет ведется на несколько часов или суток вперед до тех пор, пока не изменятся 
метеоусловия, например, направление и/или скорость ветра на рассматриваемой терри- 
тории. После этого в программу вводятся измененные метеоданные, и счет может быть 
продолжен на следующий отрезок времени. Метеоданные берутся с ближайшей метео- 
станции или из статистики для рассматриваемого сезона и территории. 

Принятие решений в кризисных ситуациях, связанных с выбросами в атмосферу 
радиоактивных и химически опасных веществ, требует оперативного прогнозирования 
развития окружающей обстановки. Наибольшее распространение получили так называемые 
гауссовы модели, эмпирически обобщающие точные решения уравнений адвекции — диф- 
фузии в неизменном по времени поле ветра. Они положены в основу существующих 
нормативных методик, определяющих распространение в атмосфере промышленных 
выбросов, практически во всем мире [2-4]. Эти модели содержат инженерные формулы 
для расчета характеристик загрязнения; коэффициенты, входящие в формулы, опреде- 
ляются из статистики многочисленных наблюдений и экспериментов по распростране- 
нию примесей. Следует отметить, что разброс статистических данных довольно велик; 
например, отклонения от средних значений для дисперсии облака примесей могут до- 
стигать нескольких порядков [5]. 

Гауссовы модели применяются для регионов размером 10 -— 20 км. На больших 
расстояниях существенно сказывается неоднородность метеорологических полей, поэтому 
даже для небольших расстояний лучшие результаты дают модели, использующие реаль- 
ные, а не усредненные метеополя. 

В последнее время для решения поставленной задачи все чаще начинают использо- 
вать эйлеров подход, в котором численно решаются уравнения переноса и турбулентной 
диффузии примесей. Для борьбы со схемной вязкостью и немонотонностью решения 
предложены эффективные разностные алгоритмы, позволяющие получать решения даже 
на персональных ЭВМ средней мощности [6-8]. 

Предлагаемая методика в отличие от других подходов позволяет задавать произ- 
вольное пространственное распределение примеси в начальный момент времени, учи- 
тывать мобильные источники выбросов и рассчитывать нестационарный процесс 
распространения примесей от точечных, линейных и площадных источников. 

Входные данные 

1. Топографическая карта с тематическими слоями о рельефе местности, типе 
растительности и подстилающей поверхности, о населенных пунктах и численности 
населения. 

2. Координаты ближайшей метеостанции и данные с нее: атмосферное давление, 
направление и скорость ветра, температура, относительная влажность воздуха и не- 
которые другие. 
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3. Координаты выбросов загрязнений и их тип: наземный пожар, взрыв, непрерыв- 
ный или периодический по времени источник примесей, взрыв над подстилающей по- 
верхностью и т.д. 

4. Для мобильного источника примесей задаются на карте узловые точки коорди- 
нат трассы, которые соединяются отрезками прямых линий, а также скорость движения 
транспорта между соседними узлами. 

Математическая постановка задачи переноса примесей 

Пусть с$(х,у,7.1) — концентрация $-й примеси, переносимой ветром в атмосфере в 
цилиндрической области С с поверхностью $, состоящей из боковой поверхности ци- 
линдра У, нижнего основания Х0 (при 7 = 0) и верхнего основания ХН (при #=Н) [9]. 

Если 0 = и! + у] + ий - вектор скорости частиц воздуха как функция х, у, 2, # (1), К- 
орты в направлении осей Ох, Оу, 02, соответственно), то перенос субстанции вдоль 
траекторий и турбулентная диффузия примесей в пограничном слое атмосферы [10] 
описывается уравнениями: 


де, ‚ (ис,) ‚ О(ие,) дс, _ 


01 Ох ду 02 
д’; 0*е д дс 
= + |+ У(2)—"|+Е, -— (Е, +Е,. с. › (1) 
Н ед? * бу |: Я. $ — (А +, )С. 
би 0 0% | 2) 
0х ду 02 
Здесь Сс; — концентрация 5-й примеси, 5 = 1, 2, ..., 4; 4 — число примесей; 


У? =и- и’ 5 и, у, м — компоненты вектора скорости; и’ — скорость гравитационного 


осаждения 5-й примеси в атмосфере; и и у%(2) —- коэффициенты турбулентной вязкости 
воздуха по горизонтали и вертикали; функции Ё, описывают пространственные источ- 


ники выбросов примесей в атмосферу, а К, и К., — коэффициенты сухого и влажного 


осаждения примесей. Второе уравнение задает условие непрерывности среды. 
К уравнениям (1) и (2) необходимо еще присоединить начальные и граничные 
условия. В начальный момент времени решение предполагается известным: 


с= со (х,у,2) при #=0. (3) 

В программе необходимо предусмотреть возможность задавать типовые конфигу- 

рации выбросов: типа «гриб», надземный взрыв («шар»), наземный взрыв («точка») и т.д. 
Краевые условия задаются в виде: 

дс, 


С, =0 нах при И, < 0, 


дс т 
=9,с,; на Хи —= 0 на Ун, 
02 74 

где а, — функция, характеризующая взаимодействие 5-й примеси с подстилающей 
поверхностью, И» обозначает проекцию вектора скорости ветра на внешнюю нормаль 
к поверхности У. 

Первое краевое условие означает, что если поток втекает в рассматриваемую 
область, то он приносит нулевую концентрацию примеси. Если же поток вытекает из 
расчетной области, то на этой части границы сохраняется вычисленное внутри области 
значение концентрации примеси вблизи границы. На нижней границе Хо примесь может 
частично поглощаться, а на верхней (Хн) — поток примеси равен нулю. 


= 0 на Х при И, 20, (4) 
дп 


5 
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Кроме того, будем считать, что вертикальная компонента скорости и’ равна нулю 
на нижней и верхней границах расчетной области и вычисляется из уравнения (2) как 
Г ди д» 
= -| —+— 4. 
‚\@х ду 
Общее описание алгоритма. Уравнения (1) расщепляются по физическим про- 
цессам [9]: 


— перенос по траекториям а и т оо Е 8 0; (5) 
01 дх ду 02 
2 2 
— турбулентная диффузия дс, ы: дс, а дс, а д |. дс, : (6) 
0: 0 20 02\ 02 
— трансформация примесей ее =. - (в, +. )с.. (7) 
, 576: 


Уравнения (6) допускают дальнейшее расщепление по направлениям Ох, Оу, 02, 
после чего получившиеся одномерные уравнения решаются известными численными 
методами [11]. Последние уравнения - (7) — можно трактовать как систему обыкновенных 
дифференциальных уравнений по времени [, в которой коэффициенты параметрически 
зависят от х, у, 2. Для их решения можно использовать практически любой известный 
метод второго порядка точности по времени. 

Скорость осаждения примесей. Считается, что скорость М’ гравитационного 


осаждения (седиментации) твердых частиц с радиусом до 20 мкм в спокойной атмосфере 
зависит только от их размера и вычисляется по формуле Стокса: 


а 
и, = (р, -р,), (8) 
187, 


где 4, — аэродинамический эквивалентный диаметр частицы (мкм), 77, — динамический 
коэффициент вязкости воздуха (Па:с), © — ускорение свободного падения (9,8 м/с”), = 
плотность частицы (кг/М”), р, - плотность воздуха (1,2 кг/м”). Специально поставленные 
эксперименты позволили дать формулы для широкого диапазона размеров частиц (табл. 1). 


Таблица 1 — Скорость осаждения частиц 


Радиус г, мкм Скорость осаждения у›., см/с 


и’; = Кг? ‚ = 1,19-10° см" с" (закон Стокса) 


40 <" < 600 и, = ДГ, Ёз= 8-10 с" 


Турбулентная вязкость. Физические свойства атмосферы по горизонтали считаются 
одинаковыми. По этой причине коэффициент турбулентной вязкости по х, у берется 
постоянным и равным 200 — 2000 м’/с. Вертикальная составляющая коэффициента тур- 
булентной вязкости вычисляется по модели Иорданова [9] или по формуле Извекова [12] 


У(2) = +Ут а | (9) 


Параметры имеют следующие типичные значения: \0 = 5 м/с, у„ = 0.266-10^, 
т = 0.02 м". 
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Вычисление функционалов от решения. Необходимо вычислять следующие функ- 
ционалы: 
1. Количество 1-й примеси в атмосфере А (область С’) в тоннах или мегабек- 


керелях (МБк) 
А' = [с'а6 (10) 
б 
или в разностном виде 
_ КАМАМ- - 
А = оз 2 р. АР а — 2,1). 
= т= и=1 
2. Суммарное количество всех примесей в атмосфере 
А= — РК (11) 


(м > 2 
3. Плотность загрязнения 1-и примесью подстилающеи поверхности в мг/м” или 
2 
МБк/м 


1 
Р' (ху, 1) = [е'(х, у,бГ)аг (12) 
о 


ИЛИ Е (1 + 01) т Ра () + же В 
где О = та о (1) + О (1 ке 91) р “ 
Вначале (при {# = 1) в т.0 Полагаются равными нулю, & — время начала выброса 


примеси. 
4. Общая плотность загрязнения подстилающей поверхности 


а (13) 


5. Вычисляются максимальные значения Атах и Риах, а также точки (хХтах, Утах» Ртах). 
в которых эти максимумы достигаются. Если какой-либо из максимумов превышает 
ПДК, то об этом выдается сообщение в легенду красным цветом. 

В случае выброса химически опасного или сильно действующего ядовитого вещест- 
ва вычисляются зоны смертельного поражения человека. В моменты времени, кратные (й, 
рисуются поля концентрации примесей на фоне карты Беларуси и выдается легенда. 


Должно быть предусмотрено рисование изолиний функций }(х, у) = с, в двух 


режимах: 


— относительные уровни (5 изолиний) с, = Ц ,с 4.— константы, а Стах — макси- 


тах 
мальное значение функции /(х, у). Константы 4 должны быть легко сменяемые, по 
умолчанию (= 0,05; 0,25; 0,5; 0,75; 0,95. 

— абсолютные уровни задаются шкалой концентраций или уровнем токсодозы. 

Реализация 

В рамках выполнения государственной программы МЧС 02, в ОИПИ НАН Беларуси 
выполнена разработка программного комплекса РИСК-1. Этот комплекс предназначен для 
решения задач анализа рисков и моделирования аварийных ситуаций на химически опас- 
ных объектах. 

Программный комплекс РИСК-| сделан в виде ГИС, которая дает пользова- 
телю следующие возможности: 

— загружать и отображать векторные (стандарт Е205) и растровые (ВМР, ГРС, 
ТПЕЕ и др.) карты; 

— производить наложение и привязку к карте спутниковых и аэрофотоснимков; 
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— использование внешних БД позволяет пользователю легко менять исходные 
данные, адаптируя ГИС к различным ситуациям и условиям; 

— решать задачи оценки индивидуального риска на объекте и оценки последствий 
аварийной ситуации на объекте; 

— представлять результаты решения задачи в виде иллюстрированных отчетов. 

Для работы с электронными картами используется библиотека функций, разрабо- 
танная в ОИПИ НАН Беларуси. Она позволяет выполнять стандартный для комплексов 
такого уровня набор действий с векторными картами - загрузка, отображение, вывод на 
печать, выделение объектов и слоев и прочее. 

При запуске программы происходит автоматическая погрузка последнего вариан- 
та работы: загружается карта, с которой происходила работа, устанавливается последнее 
значение масштаба и координат отображения, подгружаются данные из нескольких баз 
данных, в списке объектов выбирается объект, который был выбран при предыдущем 
окончании работы комплекса. 

Базы данных, из которых при старте загружается необходимая информация, орга- 
низованы в виде таблиц Ехсе|, что позволяет пользователю легко и удобно просматри- 
вать и модифицировать их содержимое. В БД хранится следующая информация: 

— список СДЯВ и их физико-химические свойства; 

— список промышленных объектов, на которых имеются запасы СДЯВ, и данные 
об этих объектах, необходимые для проведения расчетов; 

— усредненная метеорологическая информация, необходимая проведения расчетов 
при решении задачи оценки рисков. 

В ходе работы, выбрав объект и введя ряд данных, пользователь может решить 
две задачи: 

— оценки индивидуального риска и риска аварии на объекте; 

— оценки последствий аварийной ситуации на объекте. 

При решении задачи оценки индивидуального риска производится расчет послед- 
ствий выброса по 8 направлениям и происходит усреднение полученных результатов по 
вероятностям направлений. Результатом является круговая диаграмма, которая показы- 
вает на карте места, где риск и вероятная токсодоза превышают определенные пороги. 

По итогам решения задачи оценки индивидуального риска РИСК-1 выдает сле- 
дующие результаты (рис. 2): 

— вероятностное поле токсодозы, 

— вероятностное поле индивидуального риска. 


ЕЕ ЕЕ ИЕ 


Сы) аа] ва Ее ЕЕ ИС 


Поле токсодозы при расчете оцезмки риска 


Рисунок 2 — Вероятностное поле токсодозы при расчете оценки риска 
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При решении задачи оценки последствий аварийной ситуации на объекте рассчи- 
тывается ряд результатов, такие, как концентрация и токсодоза при прохождении пер- 
вичного и вторичного облака, общие концентрация и токсодоза, площади зон заражения, 
оценка риска и математическое ожидание числа пораженных людей. Все это выдается в 
виде результатов на карте, графиков и таблиц (рис. 3). 


р 9) 218} ГО ©] ео] 


[Поле общей тожсодозы при прохождении облака 
| 5 


Рисунок 3— Распространение токсодозы при аварии на объекте 


По итогам решения задачи оценки последствий аварийной ситуации РИСК-1 
выдает следующие результаты: 

— поле концентрации СДЯВ в результате прохода первичного облака; 

— поле токсодозы в результате прохода первичного облака; 

— поле концентрации СДЯВ в результате прохода вторичного облака; 

— поле токсодозы в результате прохода вторичного облака; 

— обобщенное поле концентрации СДЯВ; 

— обобщенное поле токсодозы; 

— график величины обобщенной концентрации СДЯВ по оси выброса; 

— график величины обобщенной токсодозы по оси выброса; 

— график величины индивидуального риска про оси выброса. 

Используя управляющие кнопки и поля ввода, пользователь может выбрать необхо- 
димый и наиболее удобный ему режим отображения результатов (рис. 4). 


= тЫ ты к ОНИ 
{Поле общей концентрации при прохождении облака | 


Рисунок 4 — Отображение результатов форме графиков и таблиц 
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Все эти результаты пользователь может как отобразить на экране монитора, так и 
поместить в иллюстрированный отчет, который автоматически генерируется комплексом 
по запросу пользователя. 

В дальнейшем планируется расширить возможности комплекса РИСК-1 по 
поддержке других форматов карт, решаемых задач и представляемых результатов. 


Выводы 


Разработанные в ОИПИ НАН Беларуси инструментальные средства оперативного 
анализа местности для построения систем принятия решений при чрезвычайных си- 
туациях прошли опытную эксплуатацию в качестве импортозамещающего продукта и 
показали свою эффективность при решении актуальних народно-хозяйственных задач. 
Созданные методы и технологии использованы при нахождении оптимального маршрута 
движения транспортного средства с учетом ограничений на проходимость и присутствие 
зон поражения и моделирования переноса радиоактивных и загрязняющих частиц. 
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Гнструментальн! засоби оперативного аналзу мекцевост! для побудови систем прийняття ришень 
при надзвичайних ситуащях 

У стати представлен! розроблен! 1нструментальн! засоби оперативного аналзу мсцевост! на основ! цифрових 
карт та знамю для розв’язання задач 1з моделювання, анал!зу та прийняття рипень при надзвичайних ситуащях. 
Тнструментальн! засоби шдтримують рипення таких задач, як аналз властивостей мецевость, розрахунок 
найкоротших маршрут, враховуючи задан! обмеження, пошук найближчих об’ектв В заданими характери- 
стиками, визначення просторових вдношень об’екив мисцевость розрахунок зон затоплень, моделювання 
перенесення радюактивних частинок, моделювання поширювання пожеж та ряд 1н. Наводяться приклади 
використання. 


А.М. КгуисйКотъ, Г.М. 5оБо, 5.Г. АМатеуко, О.А. СаутПуик, С.Р. Аратт 

Тод Гог Орегайуе Теггат Апа!у$15 Гог Сепегайоп 07 Оес$юп-МаКшо Зу$етл$ а Етегоепсу ЗНиаНоп$ 
п #55 рарег ве деуеореа {0015 Юг Фе орегануе {еггат апа1уз15 оп Фе Баз1$ оЁ ЧеНа!| саг4$ ап тетое зеп$те 
ппаоез Кг фе 4ес151юп оЁ ргоМепл$ оп тодейпе, апа[уз15 ап Чес1юп-таКше аё етегоепсу зИиаайопз аге 
Чезспеа. ТВе 1001$ шсаде фе 4ес1з1оп оЁ зисЬ ргоепл$ аз Фе апа!уз1 оРфеташт сБагасет$Нс$, са1сШайоп оЁ 
Фе зВоцезЕ гошез ЧаКше шю ассоцие Фе зе{ гезблсНопз, зеагсВ оР Фе пеагезЕ об]есё5 УИ Фе тазег даа, 
деНиюоп о# зраНа| танопз$ оЁ1еташ об]ес6, сасайоп оЁ 7топез оЁ Ноои®, тодеНпе оЁ сатуше оуег оЁ 
таЧоасйуе рагисез, подейпе оЁ Ч1за1Бийоп о Ягез ап4 а питабег оР офегз. 'ТБе ехатар/ез оЁргасйса] гези$ аге 
ЗВо\уп. 


Статья поступила в редакцию 22.06.2011. 
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